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Abstract. The entropy of the two-dimensional t-J model is investigated using its 12th order 
high temperature series. A direct Fade extrapolation of the entropy series doesn’t converge 
well for temperatures below T ~ J. The series coefficients are exact polynomials so the series 
convergence can be improved by modifying the series that is extrapolated. By subtracting 
a scaled version of the series for the entropy of the Heisenberg antiferromagnet from the t-J 
entropy series the low temperature convergence is greatly improved. Using this technique results 
are obtained for the full range of electron densities and temperatures. The electron density is an 
adjustable parameter in the series coefficients allowing the density dependence of the entropy 
and the density derivative at fixed temperature dS/dn\T to be determined accurately. The 
density derivative depends strongly on temperature, unlike noninteracting models. The density 
derivative is also an approximation to the experimentally measured thermopower. 


1. Introduction 

The two-dimensional t-J model has been studied for many years as a model for the copper oxide 
planes found in high temperature superconductors. Despite considerable effort the basic 
thermodynamic properties of this model are still not well understood. The calculation reported 
here investigates the entropy by means of its high temperature series calculated to 12th order 
in inverse temperature (5. In principle, knowing the full temperature dependence of the entropy 
and the value of the free energy at one temperature is sufficient to fully determine the free energy 
and thus all the thermodynamic properties of the model. This more general calculation is left 
for a future publication. In this paper our attention is restricted to the entropy per site S and 
its density derivative at fixed temperature dS/dn\T- 
The t-J model Hamiltonian is given by 

H = -tY^ (cLcjA + Jj'^Si- Sj, (1) 

^{ij) (U> 

along with the constraint of no double occupancy. The high temperature series for the Helmholtz 
free energy per site is calculated by a linked cluster expansion, taking the form 

N 

PF = Y,Mf3Jy, ( 2 ) 
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where the fi{n,t/J) are exact polynomials in terms of the electron density n and the ratio of the 
model coupling constants t/J. For the calculation reported here = 12 and we fix t/J = 2.5. 


The series for the entropy per site is also calculated in terms of the fi as 


N 

5 = A:b 
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and the energy per site is 

N 

l3E = Y,ifiWJr- 
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2. Details of the Calculation 

To reach low temperatures, the series expansion for the entropy needs to be extrapolated by Fade 
approximants. With the available 12th order series the direct extrapolation of the t-J entropy 
only converges for T > J. A low temperature scale remains unresolved in this calculation. To 
improve the convergence we can try to extract this low temperature scale by forming the series 
for 

AS{n,T) = Stj{n,T) - n^pSAFiJ^T) - a*STB{n,T), (5) 

where Saf is the entropy of the Heisenberg antiferromagnet (found by setting n = 1 in the 
t-J entropy), Stb is the tight-binding model entropy and the starred parameters are adjusted 
to produce a function AS with a simple temperature dependence that extrapolates well to low 
temperatures. A least squares fit of the [6/6] Fade to a fourth order polynomial is calculated 
to extend AS to T = 0. For n > 0.55 a reasonable AS can be found by setting n*^p = re, 
a* = 0 and adjusting J*. To extract Stj we need to know SAFiT) which is determined by 




Figure 1. a) Heisenberg model entropy vs. temperature, b) An example of the fitting procedure 
described in the text for the t-J model entropy at re = 0 . 88 . The t-J entropy is the sum 
Stj = nSAFiJ*) + AS, with J* = 0 . 638173 J . 

using a combination of high precision Monte Carlo data and series expansion resnlts 

for Eaf{T)- The entropy is found from the energy by using the thermodynamic relations 
E = F + TS and Faf{T) = Eqaf — SAF{T')dT'. The integal is evaluated by the trapezoidal 
rule with step size AT = 10“^J and the temperature dependence of the energy is given by 
Eaf{T) = Eaf{T) — Eqaf with ground state energy Eqaf = — 0 . 669437 J [^. Using these 






expressions the temperature dependence of the entropy Saf{T) is determined iteratively by 


sMi^T) = I (6) 

where j runs from 1 to 500,000 to cover the range 0 < T < 5J and Eaf{0) = 5'af( 0) = 0. 
Fig(la) shows the result for the temperature dependence of Saf{T). Using the relation 
above between Faf(T), Saf{T) and the value Faf{5J) = —3.5043614J determined from 
the series expansion the ground state energy found by integrating Saf{T) given by EqQ is 
^OAF ~ —0.66947J and F^^p — Fqaf = —3.3 x 10 “® J, comparable to the error estimates in the 
Monte Carlo data (^[^. 






Figure 2. A selection of data plotted as S/T vs. n. a) 0 < T/J < 0.45 in steps of 0.05. The 
dashed line in the inset is the tight-binding model 7 . b) 0.5 <T/J< 1.5 in steps of 0.1. c) 
1.6 < T/J < 2.6 in steps of 0.1. d) 2.8 < T/J <5 in steps of 0.2. 

The temperature dependence of the Heisenberg entropy Saf{T) can be used as a known 
function to reverse the AS calculation and extract Stj{n,T). Fig(lb) gives an example of this 
procedure for n = 0.88. Two temperature scales emerge from this analysis: a high temperature 
scale T ^ t and a low temperature scale T ^ J*. For n < 0.55 the t-J entropy crosses over 
to the tight-binding model entropy and all the starred parameters in Eq([^ are required. A 
complete discussion of this more complicated htting procedure at low n is deferred to a future 
publication, though here we note to order n the series coefficients for the t-J model free energy 
































are exactly the same as the tight-binding model free energy so as n —>■ 0 we have Stj Stb at 

all temperatures. 

3. Results 

The extrapolation in temperature to find Stj{n,T) is done separately for each n. For 0.10 < 
n < 0.92 the spacing of densities is An = 0.01, while tighter spacing is used for the highest and 
lowest density ranges. Data were accumulated for 226 densities, with 5, 000 temperature values 
for each density with a uniform spacing of AT = O.OOIJ. A selection of the entropy data is 
shown in Fig(2), with Stj/T plotted as a function of n for a range of temperatures. 

From the entropy data an interesting quantity to calculate is the density derivative at fixed 
temperature dS/dn\T- This derivative has been investigated in earlier calculations [^, but 
here we have greatly improved density resolution allowing a more detailed investigation than 
before. The entropy initially grows upon doping away from half hlling and at high temperatures 
the configurational entropy goes through a maximum at n = 2/3. For noninteracting systems 
the configurational entropy maximum remains true for all temperatures. For example, the 
maximum entropy for the tight-binding model is found at n = 1 for all temperatures. Even 
for hard core bosons the entropy maximum remains fixed at n = 0.5 for all temperatures. The 
entropy maximum for the t-J model is strongly temperature dependent as shown in Fig(3). 
As the temperature decreases from T ~ 5JtoT ~ J the entropy maximum moves up to 


T/J 



Figure 3. Red line: dS/dn = 0 on a temperature vs. density plot. Green line: configurational 
entropy maximum at n = 2/3, independent of temperature. Blue line: approximate temperature 
scale T = t for the deviation of dS/dn = 0 from the configurational entropy maximum. Magenta 
line: T = 0.2J, approximately room temperature where dS/dn = 0 at n = 0.796. 


n = 0.84. In this temperature range the interactions in the t-J model are decreasing the 

















Figure 4. Plot of dfi/dT vs. temperature for the density range 0.7 < n < 0.95. The spacing of 
the densities is 0.01. 


entropy at higher densities than would be expected from the configurational entropy. This is an 
indication of correlations developing in this density and temperature range. The temperatures 
are too high for these to be antiferromagnetic correlations and the density and temperature 
ranges don’t match the pseudogap found in high temperature superconductors. For T < J the 
entropy maximum reverses and moves back to lower densities, with the temperature dependence 
decreasing sharply down to n = 0.8 and a low temperature tail below T = 0.2J extending down 
to re = 0.62. Shastry has shown that dS/dn is an approximation for the thermopower. In 
experiments on cuprate superconductors the thermopower is observed to change sign as a 
function of doping. In particular the thermopower is zero at room temperature for 23% doping 
in many cuprate superconductors [^[^. For T = 0.2 J (corresponding to room temperature) the 
t-J model entropy maximum is at 20.4% doping. 

Further information on the correlations developing at lower temperatures in the t-J model 
can be found by considering the full temperature and density dependence of dS/dn. By the 
equality of the mixed partial derivatives of the free energy we also have 

&f 

where /U is the chemical potential. The relation in EqQ provides another means to interpret 
the derivative. Fig(4) shows d^/dT as a function of temperature for the range of densities 
where dS/dn = 0. The low temperature positive peak in dfi/dT for re > 0.8 is due to 
the attractive antiferromagnetic interaction between oppositely oriented spins decreasing the 
chemical potential at low temperatures. It costs less energy to add an electron to the system 
when the temperature falls below the effective spin interaction energy. For re < 0.8 in Fig(4) 


dn 


(7) 











the opposite happens: a fairly broad negative peak in dfi/dT develops. Below this peak it costs 
more energy to add an electron to the system than for temperatures above the peak. For the 
full range of densities shown in Fig(4) the high temperature limit of dfi/dT is positive. 

4. Conclusions 

The strong temperature dependence of dS/dn = 0 shows the interactions in the t-J model are 
producing competing correlations. For n > 0.84 antiferromagnetic fluctuations are dominant 
at low temperatures. For n < 0.8 the dominant fluctuations at low temperature are d-wave 
pair fluctuations as shown in Ref. For densities 0.8 < n < 0.84 there is a crossover between 
the two dominant fluctuations. The different fluctuations also give very different temperature 
dependences to dfi/dT. The antiferromagnetic fluctuations reduce the energy needed to add an 
electron at low temperatures while the d-wave pair fluctuations increase the energy needed to 
add an electron at low temperatures. 
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